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We give a detailed description of the differences between the factorization and results derived 
from SCET and QCDF for decays B — > M1M2. This serves as a reply to the comment about our 
work B — > M1M2: Factorization, charming penguins, strong phases, and polarization Q made by 
the authors in 0. We disagree with their criticisms. 
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In Q we derived a factorization formula for exclusive B 
decays to two light mesons using the soft collinear effec- 
tive theory (SCET) 0. Recently, Beneke, Buchalla, Neu- 
bert and Sachrajda posted a comment about our work 0, 
and compared it with their QCDF (QCD factorization) 
approach |4| . In thispaper we compare results and reply 
to their comments 

For easy reference, we summarize a few points made 
in Ref. [lj that disagree with Ref. We found that: 
i) a proper separation of scales Q 2 ^> E n A ^> A 2 in the 
factorization theorem are different in SCET and QCDF, 
and can be regarded as a formal disagreement if desired 
(Q = 771;,, E-k); ii) only a subset of a s (rrib) corrections are 
currently known, so results for these corrections are for- 
mally incomplete; iii) certain amplitudes are sensitive to 
the treatment of m c , with parametrically large ~ v con- 
tributions from cc in the NRQCD region where v is the 
velocity power counting parameter; iv) current B — > tttt 
data analyzed at LO in SCET supports values for the 
parameters C S7r ~ Cj^; m disagreement with numeri- 
cal inputs adopted in QCDF. We also emphasized that 
A/m,b power corrections need to be of natural size to sup- 
port model independent phenomenology and verify that 
the expansion converges, concepts which are sometimes 
relaxed in QCDF phenomenological analyses. If power 
corrections change the LO values of £ B7r and v sub- 
stantially then this would indicate that the power expan- 
sion is not converging and a different expansion would be 
needed if model independent results are desired. 

We take this opportunity to also comment on results 
from Ref. where we found agreement with points made 
in Ref. Q|. The original starting idea is the same, that 
factorization theorems for these decays should be derived 
by making a systematic expansion of QCD in AQcn/mb, 
where the terms in this expansion are model indepen- 
dent and unique. Earlier discussion of QCD based fac- 
torization methods for nonleptonic decays can be found 
in 0, |(|. We agree on the scaling in A/ mj, for the 
LO amplitudes. There is agreement that input on non- 
perturbative functions can be obtained from B — > tt form 
factors and LO light-cone meson distribution functions 
<j>Tr(x) and 0s(fc + ). We also agree that at LO factor- 
ization occurs for amplitudes from light quark penguin 
loops, as well as tree, and color-suppressed diagrams. Fi- 



nally, the set of the one-loop hard corrections computed 
in Ref. determine the Wilson coefficients of the 
operators 0- 

It is also worth emphasizing that the scope of our two 
works was different. In Refs. 4, 8] factorization theorems 
were proposed based on the study of the IR singulari- 
ties of lowest order diagrams in perturbation theory. In- 
put parameters are taken from QCD sum rules. Certain 
power suppressed contributions were also included, al- 
though the factorization was not extended to this order; 
as a result, some of these corrections are IR divergent, 
and cutoffs were used for numerical estimates. In this 
way the authors of Refs. 0,0 were able to make predic- 
tions for many modes, which however depend on model 
dependent input. In contrast, in ,^ we used operators in 
SCET to separate the long and short distance physics to 
all orders in a s . The decay amplitudes factor, with the 
long distance physics given by a few universal hadronic 
parameters. Predictive power was shown to be retained 
even when a s (\/EA) effects are summed to all orders. 
We then used data to determine these LO hadronic pa- 
rameters and obtained a prediction for the heavy to light 
form factor /+(0). The factorization theorem also gives 
a model independent determination of the weak phase 7 
(or a) using Br(B — ► 7r°7r°) as input, but not C^o^o |9j- 

We organize the remainder of this paper by the four 
points raised in |2( which we disagree with: 1) that the 
SCET result is formally equivalent to the QCDF re- 
sult 0, in all respects, 2) that there is little bene- 
fit to avoiding the perturbative expansion at the scale 
l 1 = \/Aqcd TTtbt 3) that operators containing a charm 
quark pair are perturbatively calculable, with corrections 
suppressed by A/m& regardless of how the scale 777 c is 
treated, and 4) that our phenomenological analysis of re- 
cent B — > 7T7T data, which disfavors certain QCDF input 
parameters, is flawed because it omits "known" pertur- 
bative and power suppressed contributions. A section is 
devoted to each of these topics. 



I. FORMAL COMPARISON 

In SCET the separation of scales Q 2 > E^A > A 2 
can be achieved by matching QCD onto a theory called 
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Separating the E n A 3> A 2 scales is more complicated. 
For Cj^(z) we found 



FIG. 1: Examples of QCD diagrams. In Eq. Q a) and b) 
both contribute to the Cf n { z ) term in the factorization theo- 
rem, since the _L-gluon in b) decouples from the upward go- 
ing pion and can be Fierzed into the other quark bilinear. In 
Eq. graph a) contributes to the form factor term, and b) 
contributes to the hard-scattering term. 

SCETi to integrate out Q 2 E3> 

and then matching 
SCETi onto a final theory SCETn to integrate out the 
scale E„A E3- Performing the first step for A„+„- = 
A(B — > 7T + 7r - ) at LO in the power counting we found 



A^*- =N^f v jdudzTv{u,z)Cf*{z)4r(u) (1) 
+ / Jdu T K (uW («) | + A { f) AH , 



where N = Gprn? B /y/2, the Tj's capture hard a s {mb) 
contributions, and Cj^j C Btc depend on the \fE^A and 
A scales. The analogous expression for the B — > ir form 
factor is 



/+(0) = rW + ldzC { +\z) Cf"(z) , 



(2) 



from which we observed that both observables depend 
on the same universal Cj*( z ) and ( B ™ ■ Currently Tic(u) 
is known at 0(a s (nib)) 3, but T\j(u) is not, hence our 
statement that the calculation of the hard aJrnb) cor- 
rections are incomplete. In the notation in [lj the one- 
loop matching for the Wilson coefficients b/\u,z) are 
missing. We do not believe that these facts are disputed 
in . The amplitude A™ denotes long-distance cc con- 
tributions which we take up in a separate section. Short 
distance cc can contribute to Tij and T\q. 

In Ref. t 4] the QCDF factorization formula was 



A^-=Nf + (0)f n JduT I (u)(lf(u) (3) 
+N SlfB [ du dx dk+T u (u, x, fc+)^(u)^ {x)4> B (fc+) . 



Both of the scales nib and v EA are treated perturba- 
tively in T\\, so this result does not formally distinguish 
between these scales. This makes it impossible to work 
to all orders in a s at the u, ~ \J EA scale. It is also not 
possible to sum logarithms between nib and \[E^A with- 
out further factorization of T\\. Fig. Ogives an example 
of how individual diagrams are treated differently in 
and as explained in the caption. The result in Eq.QJ 
separates out the hard contributions ~ Q 2 regardless of 
the sensitivity to smaller scales y/E^A and A. 



UJb jdk + jdx J(z,x,k + )<p n (x)^(k + ) , (4) 



where the jet function J starts at 0(a s (\/ E^A)) and 
captures all corrections at this scale. (Multiplicative 
a s {nib) corrections in fg can be moved into T\j if de- 
sired.) The jet function J is now known at one-loop 
order E2- As first discussed in 0, using SCET one can 
write Tn(u, x, k + ) = Jdz[Tu(u, z) — Cj {z)]J{z, x, k + ); 
a point mentioned in Q on which we agree. A full dis- 
entangling of the scales E„A and A 2 in ( B7r is still being 
debated [TEH El El El. In SCETi the diagrams that 
define £ Bir involve the exchange of at least one hard- 
collinear gluon E3j J us t as they do for Cj W i leading to 
the expectation that this parameter should also start at 
0(a s (y/E~A)), and C, B7r ~ C, B * . Th e parameter count- 
ing in QCDF assumes an ct s (v 'E^A)) only for the hard 
spectator contributions, and in our notation they have 
( B7T <C £ S?r - The theoretical issue that blurs the an- 
swer to this question is that naively performing similar 
steps for ( B7V to those giving Eq. (0J give divergent con- 
volution integrals, indicating that this parameter is "non- 
factorizable" HEBE1E3]- A complete understanding of 
the separation of the ^JE^A and A scales is necessary to 
sum all logs below the hard collinear scale. For the con- 
tribution in Eq. Q these logarithms have been resumed 
in Ref. 1181 , and the running for the other part was stud- 
ied in 14]. It has been argued Ell that the a s {^/E^A) is 
absent in C Sti \ however this relies on the conjecture that 
diagrams containing a soft-collinear messenger mode E3 
in the theory below the scale y/E^A cancels all endpoint 
singularities. 



II. PERTURBATIONS THEORY IN a s (s/EA) 

In we showed that an expansion in a s {\fB^K) is 
not required to obtain predictions from factorization in 
nonleptonic decays. In other words Eq. ^ has nontrivial 
implications even without using a perturbative expansion 
for J in Eq. Q). Since y/E^A — 1.0-1.6 GeV, we believe 
that it is useful to consider what predictions can be made 
by avoiding this expansion. 

In our opinion the advantage of using Eq. Q with an 
expansion of J in a s (VKA) is that it reduces ( BM {z), 
which depends on the choice of M, to two more universal 
nonperturbative functions 4>~^(k + ) and 4>m(x). Thus, the 
number of unknowns decreases in global analyses involv- 
ing different light mesons M, and processes like 7*7 — > M 
and B — > -feP |l9j can be used to provide additional non- 
perturbative information. This is true both at tree level 
in a s (m,b) and at one-loop in a s (nib). 

There were two main critiques raised in 0. First, it 
was argued that there is no point in working to all orders 
in a s (\/E n A) since perturbation theory in this quantity 
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is reasonably well behaved, and second that if one works 
to all orders in a s (y/ 'E^A) then the factorization theorem 
looses all predictive power beyond tree- level in a s (m&). 

Regarding the first point, it is always better to use 
less theoretical assumptions if predictions at the level 
of precision currently achievable actually do not rely on 
making this expansion. The perturbative expansion at 
the scale WrnbAQCD may in fact converge (recent evi- 
dence has been given by the one-loop calculation in ^1 
and also relations between B — ► Drj and B — > Drf de- 
cays however if a prediction is independent of this 
expansion then there is no need to rely on it. Regarding 
the second point, it is true that beyond LO in a s (mb) 
the functional dependence of (f n (z) is required, instead 
of just the number £j j7r = Jdz(f n (z). However, the same 
function determines the B — > n form factors. Note that 
when a s corrections are included in Eq. J2J the functional 
form of the B meson wave function </>s(fc + ) is required. 
So in cither case one has moments of one unknown func- 
tion. We believe that the most important advantage of 
the a s (V m& A) expansion is the universality of the non- 
perturbative functions mentioned above, rather than the 
change in how a s (mb) corrections are included. 



III. CHARM LOOPS AND A c <, 

The size of charm loop contributions to B — > tttt 
are important. Certain charm loop contributions are 
from hard (~ nib) momenta and there is broad agree- 
ment 0, 0, El that these effects can be computed 
in perturbation theory. In Eq. they enter in both 
T\q and T±j. The point being debated is the parametric 
scaling of non-perturbative contributions from penguin 
charm quark loops (so-called charming penguins |22jp. 
denoted by A£g in Eq. JTJ. 

For the parametric scaling two useful limits are 



m r A 
i: — <1, — <1 

m b m c 

ii): 2^-0(1), A<1. 
m b m c 



(5) 



In this corresponds to the limits rrib —> oo with m c 
fixed and m,b,m c — ► oo with m c /mb fixed respectively, 
however we believe the description in Eq. J^J) makes as- 
pects of the expansion more clear. For example, the 
charm quark power counting will not be identical to that 
for light quarks unless m c ~ A, which is not realized in 
nature. 

In Ref. we focused on nonperturbative contribu- 
tions from cc in the NRQCD region and found that these 
contributions are only suppressed by v. Note that these 
cc's can still have a total energy ~ nif, as long as their 
relative velocity ~ v is small. 1 Here v is a place holder 



for a nonperturbative matrix element which is paramet- 
rically of this size. In charmonium m c v ~ 800 MeV and 
m c v 2 ~ 400 MeV , so if we identify one of these scales 
with Aqcd the v suppression becomes either a Aqcd /m c 
or a Aqcd/^c suppression. We do not think that 
the case m c v 2 ^> Aqcd is physically relevant for charm 
quarks in QCD. In either case if we expand in A/m c 
these nonperturbative contributions do not enter the LO 
B — > 7T7T factorization theorem. However, in practice 
v ~ 0.5 so these corrections are numerically large com- 
pared to A/m b ~ 0.1 power corrections, and can spoil 
the power expansion. 

In Refs. 0, H 13 it was argued that the NRQCD cc 
region does not require special treatment due to quark- 
hadron duality, with smearing from the q 2 = xm 2 of the 
gluon which the charm annihilate into (0 < q 2 < ml). 
Using duality in this sense requires an inclusive hadronic 
final state, to make it possible for there to be a cancella- 
tion of infrared divergences between the virtual and real 
diagrams to all orders in a s . For exclusive decays like 
B — * 7T7T one must instead prove a factorization theorem 
to separate hard and infrared contributions. In these 
proofs one must consider the contributions from all pos- 
sible momentum regions. 

The arguments in [2( assume that the size of the non- 
perturbative cc terms can be estimated based on regions 
of phase space in q 2 , taking limits based on a factorization 
formula analogous to To the best of our knowledge 
it has never been proven that the NRQCD cc contribu- 
tions factor in this way. Intuitively we expect that they 
will not. The gluons whose wavelength is ~ A do not 
decouple from the charm pair which are created and an- 
nihilated in the octet state. Since the cc production and 
annihilation occur over a distance scale ~ Aq CD , the soft 
gluons radiated from energetic quarks produced from the 
annihilation may not cancel. This can lead to two types 
of Wilson lines, 5^i[n--4 S oft] and ls[n-A so ft] in the soft B- 
matrix element. In this case the amplitude will involve a 
new nonperturbative function which has a strong phase 
from the mechanism found in 25] , since the soft function 
carries information about the final state through n and 
h. 

Even in the absence of a proof of factorization for , 
it should however still be possible to determine its para- 
metric dependence on m c /mb, v, and Aqcd/to{, using 
operators in effective field theories. We find 



ATTTT 
ATTTT 



,(2m c )/( 



2m c 



(G) 



in agreement with Eq.® disagrees with the result 
in |2| since there is no A in the numerator besides that 



1 Other charm modes besides the ones considered here could also 



contribute to A££ . The consideration of the charm modes above 
is sufficient to demonstrate the scaling of this nonperturbative 
contribution. 
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FIG. 2: An update of Fig. 5 in Q giving model independent 
results for , (f n , and the B -> tv form factor f+(q 2 = 0) 
as a function of 7 and \V n b\- The shaded bands show the 1-er 
errors propagated from the B — > mr data. 

hidden in v. Physically f(2m c /mb) encodes the restric- 
tion of the charm quarks to be produced with small rela- 
tive velocity (rather than for example back-to-back with 
energies ~ m^/2). The factor of v gives the remaining 
suppression for the charm quarks to be non-perturbative 
in the NRQCD region. Together these include all "phase- 
space suppression" factors, which claims were missed 
in . A derivation of Eq. © is given in the Appendix. 



IV. 



NUMERICAL VALUES OF Cj% C, BlT 



In Q the parameters £ and £j were extracted from a 
LO SCET analysis of the B -> tttt data. New data [H 
was presented at ICHEP 2004, and we give here an up- 
date of the analysis of these parameters from [l[. We 
compare our results with the most recent QCDF analy- 
sis in and address the criticism in 0|. 

Assuming isospin symmetry and neglecting the elec- 
troweak penguins, the B — ► 7T7r amplitudes can be writ- 
ten as 



A(B°^tt\°) = X^T n -X^P, 
V2A(B- -> 7r°7r-) = \i d) T. 



(7) 



This contains five independent hadronic parameters 
which can be extracted from an isospin analysis. Using 
the world averages of the data [2(| and setting 7 = 64° 
one finds the results in the last column of Table [I] and 



, _ r0.15± 0.02(1), 

I 7i I / 7r - \ .18± 0.02(H)' 

P/N n = (-0.024 ±0.007) + (0.021 ±0.007)i 



(8) 



Note that our power counting for A™ in Eq. © 
gives complex values of a similar size, |P|/AT W ~ 
a s (2m c )vT c /N 7T ~ 0.03, where we took /(2m c /m(,) ~ 1. 

Predictions in QCDF involve an expansion at the inter- 
mediate scale and use sum-rule calculations for </>s(fc + ) 
since this hadronic function is not known from data. Our 
approach was to instead fit the parameter £j j7r to the non- 
leptonic data and using (x^ 1 )^ = 3.0 which falls within 



the rang e, 3.2 ± 0.4, preferred by fits to the 7*7 — > 7r° 
data[27|. We also assumed that a s (mb) corrections will 
be of a similar size to neglected power corrections. Nu- 
merical justification for this is discussed below. 

Factorization formulas like Eq. (JTJ express the ampli- 
tudes T and T c at leading order in l/m& in terms of the 
nonperturbative parameters <^ B7r and £j*(x). Using the 
leading order SCET relations from p] we find 



-3 







I leading 
7=64° 



These can be used to predict the B — > tt form factor at 
q 2 = as /+(0) = ( ByT + (f*. In Fig. |3 we show how 
results for our extraction of the £'s an d /+(0) depend 
on the input value of 7 (with normalization taken using 
the value of \V u b\ preferred by current inclusive fits 28]). 
Note that the smaller value of /+(0) favored from our 
analysis would increase the value of \V u b\ from exclusive 
decays, perhaps bringing it in line with the inclusive anal- 
yses. 

Working at tree level in the jet function, /+(0) = ( B7r + 
77 and the (j parameter is given by 



c 



Nr. 



(10) 



With the input parameters adopted in QCDF (J, |a| at 
leading order in a s and A/m& one can find values for 
the SCET parameters. In the default scenario of Q we 
find: C S?r = 0.26, (f 71 = 0.02, while in their scenario 2: 
C B7T = 0.20, Cj* = 0.05, which are quite different from 
Eq. @. There are two possible explanations for this 
disagreement: i) higher order perturbative and power 
corrections are important; ii) some of the hadronic in- 
put parameters used in QCDF are not supported by the 
data. The authors of 0] take the first point of view and 
argue that there are large known perturbative and power 
corrections to the leading order result. 

To compare with QCDF we have calculated the 
hadronic parameters T, T c using the analysis in 0, . 
The results are shown in Table [I] for two sets of input 
parameters from [jj, their default scenario and their S2 
scenario. We have organized the terms according to the 
expansion advocated in Ref. with a s (V EA) terms in- 
cluded at LO. The first line in the table shows the LO 
terms with a s (m b ) corrections in square brackets, and the 
second line shows A/mb corrections. We have dropped 
other A/rnb terms which contribute < 6 x lO -3 ^ in 
the amplitudes T, T c for default inputs. In contrast to 
Ref. we evaluate all full theory Cj's at fx = nib, since 
below this scale the running does not follow from the 
usual anomalous dimensions of the clectroweak Hamilto- 
nian. Note that J dx<j)^{x)/x — 3.3 in the default sce- 
nario and J dx(/> 7r (x)/x = 4.2 in the S2 scenario. We give 
the corresponding LO SCET results for these two cases 
in the table. 
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QCDF default scenario 


QCDF scenario 2 


SCETlo, (x- 1 )^ = {3.3,4.2} 


data 


T 


0.354 - [0.037 + 0.019i + . . .] 


0.358 - [0.033 + 0.017i + . . .] 


1.1% B7T + {2.17, 2.44}Cf n 


|T| = 0.30 + 0.02 




+(0.006X ff ) 


+(0.016X H ) 






T c 


0.272 + [0.008 + 0.0i + ...] 


0.234 + [0.008 + O.Oi + . . .] 


0.99^ + {0.79, 0.74}C| n 


T c = 0.16 + 0.02 




-(0.016 + 0.001^+...) 


-(0.018 + 0.003X H + ...) 







TABLE I: Numerical results for the B — » irn amplitudes. The first two columns are QCDF results organized so terms with 
ct a (\fWR) are included at lowest order, and using two scenarios for the nonperturbative input parameters from ||. The numbers 
in the square brackets are the 0(a a (mi,)) corrections, and the numbers in round brackets are so-called "chirally enhanced" 
A/mb corrections. Unknown a s (wi) and A/mb corrections that contribute at the same order are indicated by ellipses. The 
third column is the LO SCET result in terms of nonperturbative parameters. The last column is the experimental data 
from an isospin analysis using 7 — 64° and \V u b\ = 3.9 x 10~ 3 . The results for all amplitudes have a common prefactor 
N n = Gfm|/ r /\/2 removed. 



From Table [I] the a s (mt) perturbative corrections 
amount to a < 10% shift in the leading order results for 
T,T C . The "known" non-perturbative corrections with no 
Xh factor are < 10%. Non-perturbative corrections pro- 
portional to divergent convolutions in QCDF are < 10% 
for the canonical choice Xh = 2.4 from Ref. [8J. All of 
these are of a similar size to the 10 — 20% corrections 
we expect from other unknown power corrections. Thus, 
with the expansion in [1| the model parameters used in 
QCDF support the claim that the leading order extrac- 
tion of C S7r and Cj 71 from T,T C is good to - 20%. The 
reason Ref. found that power corrections change * 
substantially is that their input parameters give a LO re- 
sult that is small, numerically of a similar size to a typical 
power correction. 

In conclusion, all the perturbative and power correc- 
tions which are truly known give rise to small shifts in our 
LO analysis. The criticisms of Ref. about our work are 
based on guesses about the size of hadronic parameters 
and the size of power corrections and thus in our opin- 
ion not reliable. Further work is required to get a better 
understanding of power corrections for non-leptonic B 
decays. Some interesting work in this direction has been 
done recently in Ref. [24| . 



V. CONCLUSIONS 

In this letter we discussed the differences between the 
QCDF approach and SCET approach to factorization in 
nonleptonic B — > tttt decays, expanding on the points al- 
ready made in pj and addressing the criticism in Q . We 
also commented on recent SCET work related to these 
points, which appeared after the publication of llj. 

We addressed the main points made in Ref. 0- For- 
mally, SCET tells us that either there are missing a s (mb) 
corrections in 0, @ (for ( B7 < ~ (f"), or the QCDF 
counting which treats a s {m b ) ~ a s (V EA) relies on 
( B7T ^> Cj 71 - Avoiding perturbation theory at the in- 
termediate scale Hi n t — x/Aqcd"^ might seem to intro- 
duce more nonperturbative Cj( z ) functions than expand- 
ing in a s (jjLint). However, when restricted to the subset 



of nonleptonic and semileptonic B decays into pions, this 
amounts simply to trading one unknown function for an- 
other (<pB(k+) vs. Cf^iz)). Contrary to the claims made 
in |3, we still find a complex A™ /AJ^q ~ a s (2m c )v, in- 
dicating that long distance charm penguin contractions 
can be numerically significant. Finally, we show that 
our phenomenological analysis of B — > 7T7t data, and the 
determination of the two hadronic parameters C B7T and 
Cj 71 remain correct when known perturbative and non- 
perturbative corrections are estimated as in Ref. . The 
largest corrections actually come from unknown power 
suppressed terms, but are still within our error estimate. 
We leave it to the reader to assess the relative importance 
of the agreements and disagreements. 

This work was supported in part by the DOE under 
DE-FG03-92ER40701, DOE-ER-40682-143, DEAC02- 
6CH03000, and the cooperative research agreement DF- 
FC02-94ER40818, and by a DOE Outstanding Junior In- 
vestigator award and Alfred P. Sloan Fellowship (I.S.). 

APPENDIX A: SCALING OF A™ 

In this appendix we derive the scaling in Eq. © . First 
consider limit i). Above the b-mass we have the op- 
erator (cTb)(dTc). Next we integrate out the scale mb 
and match this operator onto one with massive collincar 
charm quarks (to ensure they are moving close together, 
i.e. have invariant mass ~ 4m^) and a collinear light d 
quark, 

O prod = [Zklj'T^] [&)TT A &)] , (Al) 

Here w\ = rribU and the massive collinear charm quarks 
have n' ■ p ~ mj and p± ~ m c with the Lagrangian 
from 2jl ( we omit all Wilson lines and other fac- 
tors that are irrelevant to the power counting). These 
SCETi fields have momenta p 2 ~ to^A or m 2 . which we 
treat as the same size. The collinear expansion parame- 
ters are A = \J Ajnib and A c = m r I mt, and T = f/L will 
not contribute, but mass insertions |3Tj w iU contribute, 
such as r = ^rf±m c /{in- D) ~ A c as seen below. The 
operator therefore scales as OP rod ~ m^A 4 A^. 
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Next we integrate out the scale m c . The scale y/E^A 
is close to m c and can be integrated out at the same 
time, but the factors generated from doing this are the 
same as those from the Ogc ETl which occur in the non- 
charm contributions and so do not effect the relative scal- 
ing. Removing m c in 0P rod requires matching the charm 
fields onto NRQCD fields m \- There is an operator 
Of nn which annihilates the charm in the boosted frame, 
and at tree level comes with a l/(4m 2 ) prefactor from 
integrating out a single gluon 2 



Of od (0) 
Of an (x) 



1 

Am 2 



TT A h v ] [rfT A (a-L) X ] BK , (A2) 



x 



^T A a ± r] 



BR 



The boost matrix L [3l[ depends on whether the a is _L or 
longitudinal. The annihilation operator is _L if a s (2m c ) 



is a good expansion parameter, so we take a± in O\ lod in 
which case L *~ 1. Since [t^T ~ m 3 v 3 this reproduces 
the m 3 for Eq. IjAlf) . Power counting gives 



O prod 
ann 



m 3 v 3 



(A3) 



Label momentum conservation in Of rod implies that the 
cc total momentum has q 2 = m 2 u — Am 2 where u = 1— u. 
This gives a delta function in the Wilson coefficient, 



prod 



./ 4m 2 . \ 

»r — r) 

v mi J 



(A4) 



The 5- function is expected [2J|, and is analogous to 
Eq. (3.13) of for the factorization formula used in 
production of energetic cc state's, whose hadronization 
is governed by NRQCD. Any 5-function for Of nn just en- 
sures overall momentum conservation and can be omit- 
ted, so we count Cf nn ~ a s (2m c ). 

If we now consider the time ordered product capturing 
the NRQCD region we have 



gives Eq. JJJ (once it is integrated over the u depen- 
dence from the matrix element). Note that the scal- 
ing of NRQCD fields together with l|A4(l account for 
the "phase space suppression" factors, and that there 
is an enhancement from the charm propagators being 
close to their mass shell. We have not been careful 
to factorize the usoft gluons from the collinear fields 
etc., so (IA5() only indicates the scaling and can not be 
used to determine the factorization or final operator in 
SCETn • Integrating over a collinear matrix clement 
that depends on u, ip(u), will induce additional depen- 
dence on 2m c /mf,, so in limit (i) we have to make some 
assumption about the form of the wavefunction (which 
we have not proven) to quote a scaling. Taking the form 
"4>{u) ~ 6?2(1 — u) with 4m 2 /m 2 ~ 0.44 unexpanded gives 
A*?/A„ +7! - ~ 1.5va s {2m c ). 

In limit (ii) we can integrate out mt and m c simul- 
taneously, and the calculation above of /(2m c /m&) will 
change. Inserting j(2m c jm^) ~ 1 in Eq. © gives 
A c5/ A Tb ~ vna s (2m c ), which is the result from Q. 
In this case there are clearly (u)soft gluons that couple, 
the cc, the soft b, and soft spectator. 

Note that we have only considered a schematic argu- 
ment for the scaling of these NRQCD cc contributions 
rather than deriving a factorization theorem, so they are 
not at the same level of rigor as Q . It is possible that a 
full consideration of the factorization for A™ might un- 
cover an additional subtelty which changes our conclu- 
sions. Obviously one can not use the scaling arguments 
to fix values for A™ so we instead used experimental 
data for the penguin, P. Note that this method of anal- 
ysis is also appropriate if the power expansion for P does 
not converge. 



Cf'od^ann ^X T [Of° d (0)0?™ (x)] 



(A5) 



Cf rod Cf nn (m- 4 «- 6 )(mgm2A 4 « 3 )(n^m c A a t; 3 ) 
{m 6 b X 6 ) Cf Iod C? nn v 



O. 



(o) 

SCETi 



Our result for the relative scaling of the B — > tttt am- 
plitudes in limit i), A™ /Aj^q, is in curly brackets, and 



2 We use BR to denote the fact that the matrix elements of these 
fields need to be boosted to the B rest frame. The v scaling is 
assigned from the CM frame, but aside from L is not affected by 
the boost since there is an invariant that can be formed for the 
resulting matrix element without kinematic factors. 
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